Physical features of Adam’s Bridge interpreted from ICESat-2 based high-resolution digital bathymetric elevation model

Adam’s Bridge is a submerged ridge connecting India and Sri Lanka, generally regarded as a chain of shoals extending for ~ 29 km from Dhanushkodi on the Indian side to Talaimannar Island of Sri Lanka. A high-resolution digital bathymetric elevation model generated using the seafloor returned photons of ICESat-2 was used to understand the intricate details of Adam’s Bridge structure. Photons emanating from ICESat-2’s green laser have the potential to detect the seafloor up to a depth of ~ 40 m; taking a cue from this potentiality, in our research, we have accrued ~ 0.2 million photons representing the depth information and generated a 10 m resolution bathymetric data for the extent of Adam’s Bridge. Visual interpretations made from this bathymetric data through 3D perspectives with multi-directional lighting effects, and also with the derived parameters like contours, slope, and volumetric analysis, enabled us to recognize the current form of Adam’s Bridge’s physical features. The results from our research confirm that, in its entirety, Adam’s Bridge is a submarine continuation of Dhanushkodi and Talaimannar Island. Throughout the crest line of Adam’s Bridge, approximately 1.5 km on either side is highly undulating within the super-shallow water with occurrences of sudden depths. There is an asymmetry of transverse slopes to the base on both sides of Adam’s Bridge, indicating dominant transgression of material energy from the waters of the Gulf of Mannar compared to the Palk Strait. The volume of Adam’s Bridge computed in our research yielded a value of ~ 1 km3; interestingly, only 0.02 percent of this volume is above the mean sea level, and in general, the same is visible in optical satellite imagery—in total ~ 99.98 percent of the Adam’s Bridge is submerged in shallow and super-shallow waters.

National Snow and Ice Data Center Distributed Active Archive Center (NSIDC DAAC) (https:// nsidc.org/ data/ icesat-2/ produ cts).A Level-2A data product, ATL03, contains height above the WGS84 ellipsoid with the mean sea surface, latitude, longitude, and time for all photons downlinked by the ATLAS instrument onboard ICESat-2 observatory 39 .The ATL03 product design aims to become a single source for all photon data and ancillary information that other higher-level ICESat-2 products need 39 .Importantly, as per the ICESat-2 ATL03'S Algorithm Theoretical Basis Document (ATBD), the essential time-variable geophysical parameters will undergo various corrections that nullify the effects of ocean tides, solid earth tides, ocean loading, and solid earth pole tides 39 ; www.nature.com/scientificreports/these corrections in the context of deriving the bathymetric information signifies the applicability of ICESat-2 photon data to perceive accurate depth information and mitigate the impact of temporal variations while considering datasets acquired on different dates.
For this research to understand the intricate details of Adam's Bridge, a high-resolution (10 m) Digital Bathymetric Elevation Model (DBEM) was generated for its extent.The source of this DBEM is a point database that contains majorly ~ 0.2 million depth values contributed by ICESat-2 photons and sounding depths retrieved from Electronic Navigation Charts (ENCs) and hydrography charts.The study area contains islands like Rameshwaram and Talaimannar; the elevation values for these islands were considered from an open-access bare-earth model.Subsequently, this point database was interpolated to generate a bathymetric surface.
Three hundred ninety-six beams from seven ground reference tracks (https:// icesat-2.gsfc.nasa.gov/ scien ce/ specs#) of ICESat-2 acquired between October 2018 and July 2023 over the extent of Adam's Bridge have qualified to abide by the prerequisite conditions needed for successful bathymetric measurements; generally, these prerequisite conditions include a preference for acquisitions during night-times and seasons of low turbid load in water 41,43,54 .ICESat-2 acquisition method yields along-track data, and solely using along-track points may not produce a continuous surface through interpolation methods; thus, additional points must compensate for the distribution.For this, depth values were retrieved from published ENCs and hydrography charts available from authorities [55][56][57] .For the extent of having land part, elevation values from the open-access bare-earth model called Forest And Buildings removed Copernicus DEM (FABDEM) 58,59 have contributed to compensate for the point distribution issue.Technical validation for the seafloor depths retrieved from the ICESat-2 water-penetrated photons was done per the recommended procedures 60 using the reference datasets at two test sites (described in the methodology section).Similarly, the final DBEM was validated using the reference data available in the form of sounding depths from the hydrographic chart generated by the Indian Naval Hydrographic Office (INHO) (https:// hydro bharat.gov.in/).
The above DBEM was used to recognize and describe the physical features of Adam's Bridge-the core objective of this article.A three-dimensional (3D) perspective view aids in the visual analysis of digital surface data and is a promising approach for ocean data exploration, analysis, and understanding of seabed topography 61 .In addition, 3D perspective views with multi-directional lighting effects enhance the impressions of relief data 62 ; thus, towards perceiving the physical features of the submerged assemblage of Adam's Bridge, multi-directional lighting effects in a 3D visualization platform available in ESRI's ArcScene software (www.esri.com) was used for interpretation.DBEM can be a source for deriving contours, slopes, and aspects; these derived data act as quantitative parameters of the seafloor 63 .Thus, the physical features of Adam's Bridge are reported in this article through the interpretations from DBEM and its derived parameters.

Results and discussion
Figure 2 shows 3D perspective views of the DBEM generated for the extent of Adam's Bridge.Specifically, Fig. 2a and b are two different perspectives of Adam's Bridge, one with the observer position from the Gulf of Mannar and the other from Palk Strait.Figures 3 and 4 show a contour map of a 2 m interval and a slope map derived from the DBEM of the study area, respectively.
One important observation from the interpretation of the DBEM is that, in its entire form, Adam's Bridge is a submerged ridge with a submarine continuation of Dhanushkodi and Talaimannar Island (refer to Figs. 2 and  3) from a rationalized depth of 8 m.From the tip of Dhanushkodi, Adam's Bridge's general trending long axis is oriented from west-northwest to east-southeast in secondary intercardinal directions, entirely perpendicular to the predominant wave approach directions of both the Gulf of Mannar and the Palk Strait.The general direction of Adam's Bridge for the first 17 km of its length is about east-south-east; later, the direction is gradually curved towards the north of the east and finally touches Talaimannar Island at its east.
Along the Adam's Bridge crest, which occupies a central position from its base on either side reckoning from sea level to 8 m (treated as a base), the transverse slope towards the Gulf of Mannar mainly varies between 0.05° and 0.25°-implying that for every 10 m horizontal distance, there is a change in slope; however, the slope towards Talaimannar side is varying between 0.2° and 0.35°.In contrast, towards the side of Palk Strait, the slope, on average, is never greater than ~ 0.2° (refer to Fig. 4).The topographic surface on either side of the Adam's Bridge's crest line for an approximate length of 1.5 km is highly undulating with varying slopes and sudden depths (refer to Figs. 3 and 4).Towards the Gulf of Mannar side, from Dhanushkodi's tip, the distance from the crest of Adam's Bridge to the -8 m contour line is 6.5 km for its initial 17 km.For the further length, the distance to the base fluctuates from 2.5 to 3.5 km.Meanwhile, the distance from the crest to the base is always greater than 7 km towards the Palk Strait side.The rationale for treating the depth to 8 m is that from this level, the seafloor of Palk Strait exhibits flatness; however, from the rationalized base of Adam's Bridge, the seafloor depth in the Gulf of Mannar stretches as a low platform, deepening fair and evenly to the south at about a rate of 2 m depth per 1 km to 36 m, after which it sinks more rapidly to great depths.
Geologically, the genesis of India and Sri Lanka are closely linked, from being a part of the ancient supercontinent of Gondwana and later with the ancient supercontinent of Pangea during the Permian period [64][65][66]  the force of plate tectonics during the Triassic period, Pangea separated into two major supercontinents, Laurasia and southern Gondwana, where India, in association with Sri Lanka, emerged as one of the huge land masses in the southern Gondwana 66,67 .India drifted northwards as an isolated island in the Thethys Seas, crashing into Laurasia about 35-55 million years ago to take up its present position 67 .From the Triassic to the Pleistocene and Holocene, sea levels have repeatedly submerged or emerged the land bridge connecting India and Sri Lanka 67,68 .
The land bridge may result from conditions associated with the sea level fluctuations since the last deglaciation associated geologic and tectonic settings.Currently, this land bridge connectivity is in a zone of active wave attacks from both sides 16 , i.e., the Gulf of Mannar and the Palk Strait.Adam's Bridge's current form might be significantly influenced by the sources of sedimentation/sand and their to and fro movements, pattern of currents, and energy to move the sedimentation [69][70][71] from the Gulf of Mannar and the Palk Strait and vice versa; these parameters may have resulted in maintaining the asymmetrical slopes on either side and aggregation of fair sand deposits on its crest and sides 25 .Generally, high waves are observed in the Gulf of Mannar, particularly during the southwest monsoon period compared to the Palk Strait, and occasionally, swells dominate in the Gulf of Mannar 72 .Previously, researchers have attempted investigations to understand the wave patterns and currents that originated in the Gulf of Mannar and Palk Strait 16,17,27,72 .The wind system over India and adjoining oceanic regions generally blows from the southwest for half of the year and from the northeast during the other half; in similar lines, the rainfall pattern and its reversal process are oriented according to the southwest and northeast monsoon seasons 73 .These two Indian monsoon seasons should have a high bearing on various characteristics of the Gulf of Mannar and the Palk Strait, including their patterns of tides and currents, salinity, and chlorophyll distribution 7,27,74,75 .Adam's Bridge, being a barrier between these two water systems, will be influenced by the transgression of the material energy 76 , resulting in maintaining its current morphological structure.Generally, the West India Coastal Current (WICC) and the Summer Monsoon Current (SMC) bring Arabian Sea waters into the Bay of Bengal during the southwest monsoon, and on the other hand, the East India Coastal Current and Winter Monsoon Current (WMC) carry the Bay of Bengal waters, partly, en route Palk Strait and Adam's Bridge, into the Arabian Sea during the northeast monsoon 77 .
At regular intervals of Adam's Bridge, sudden narrow channels with depths varying between 2 and 3 m exist (refer to Fig. 5); these narrow channels probably permit free flow or exchange of water between the Gulf of Mannar and the Palk Strait.Importantly, from the crest line of Adam's Bridge, the narrow channels are accompanied The primary inference from both these perspectives is that in its entire form, Adam's Bridge is a submerged ridge with a submarine continuation of Dhanushkodi and Talaimannar Island.The exposed sand banks are only 0.02 percent of Adam's Bridge's total volume when the base is considered at 8 m depth.These perspectives were generated using ESRI's ArcScene Ver.10.8.1 software (https:// www.esri.com/).Volumetric analysis by fixing 8 m water depth as a base for the study area yielded a volume of ~ 1 km 3 .Out of the total volume of the entire Adam's Bridge, the volume ratio towards the Gulf of Mannar and Palk Strait is 44:56.Similarly, the volume of Adam's Bridge above 0 m is 0.02 km 3 which is only 0.02 percent of the total volume; this is the same extent that is visible in the optical satellite imagery.

Conclusions
This article reports the physical features of the Adam's Bridge, an isthmus mostly submerged in shallow and supershallow waters.Approximately 0.2 million photons representing seafloor depths were aggregated from various qualified ICESat-2 ATL03 data product acquisitions for the Adam's Bridge region.Additionally, sounding depths from ENCs/charts and elevation data from the open-access bare earth model have contributed to densifying the point database.A 10 m DBEM was generated by interpolating this point database.This DBEM was used in visual analysis aided by 3D perspective views with multi-directional lighting effects.Similarly, the contour and slope map derived from this DBEM and volumetric analysis have aided in quantifying the submerged features of Adams' Bridge.The report is the first to provide intricate details about the Adam's Bridge using ICESat-2 waterpenetrated photons.Our findings can aid in accentuating the understanding of Adam's Bridge and its origin.

Datasets
Table 1 shows the datasets utilized to generate a DBEM for the extent of Adam's Bridge.The primary data resource to generate the DBEM is depth information of the seafloor retrieved using water-penetrated ICESat-2 photons.Over the extent of Adam's Bridge, seven ICESat-2 reference ground tracks are available (listed in Table 1 and illustrated in Fig. 6a).From these seven reference ground tracks, one hundred thirty-three ICESat-2 data acquisitions between October 2018 and July 2023 are available for the study area.Out of these, abiding by prerequisite conditions like preferring nighttime acquisitions and omitting the data acquired during turbid load periods (refer to Table 1), 66 tracks are found useful, and these tracks comprise 396 strong and weak beams of along-track data that were considered for this research (shown in Fig. 6b).During ICESat-2 data acquisition process over the water bodies, most photons will be reflected from the surface water.However, depending on the optical properties of the water, some of the photons will return from the water column and some from the seafloor 41,43,48,49,78,79 .Towards retrieving the depth values from the seafloor, preference should be given to those acquisitions of ICESat-2 during clear water (low turbidity) periods 41,54 .For inferring the turbidity load, a data layer titled K d (490), available from the Sentinel-3 A/B Ocean and Land Colour Instrument (OLCI) as a Level-2 series of data product service helps to characterize the transparency of water [80][81][82] .Acquisitions during high turbid load in the water were avoided by referring to the K d (490) of Sentinel-3 A/B; during this crosschecking procedure, the overlap period between ICESat-2 and Sentinel-3 A/B acquisitions is kept within + / − 24 h.Similarly, for water depth-related studies using ICESat-2, it is recommended to prefer nighttime acquisitions because the background noise caused by solar spectral radiation significantly impacts the depth detection performance of LiDAR 83,84 .Additionally, to increase the density and the distribution of the points representing the seafloor's depth in the regions with less ICESat-2 photons (also, due to the along-track pattern of ICESat-2's data collection), sounding depths available from ENCs and charts were considered and their details are given in Table 1.Similarly, to increase the point density over the extent of land (Rameshwaram and Talaimannar islands), elevation values were considered from the open-access bare-earth model called FABDEM 58,59 .Various reference datasets were used in the accuracy assessment of the depths retrieved from the ICESat-2 water-penetrated photons and the output DBEM; the details of these reference datasets are given in Table 1.

Methods
Operational methodologies to derive the seafloor depths from ICESat-2 photons are available from the earlier researchers' works 29,35,41,[47][48][49]52,[86][87][88] , which we have referred to and implemented in this research. We have onsidered only those ICESat-2 beams that have abided by the prerequisite conditions, like preferring nighttime acquisitions and omitting the data acquired during turbid load periods.All the qualified beams of data were processed to classify the return photons from the water surface, water column, land, and seafloor using Density-Based Spatial Clustering of Applications with Noise (DBSCAN) 84,[89][90][91] followed by manual correction using localized statistical algorithms 29,91,92 to eliminate outliers.Typically, the DBSCAN algorithm, by taking parameters like radius and minimum points, will classify a set of points within the said radius as signals when the density of points exceeds the pre-set threshold of minimum points 29 .The concept of retrieving the seafloor returned photons is illustrated in Fig. 7. Figure 7a shows a subset of the ICESat-2 beam acquired over the extent of Adam's Bridge, and Fig. 7b shows its corresponding 2D profile consisting of all the geolocated photons (including noise) as recorded in the ICESat-2 Level-2A ATL03 data product.Typically, these 2D profiles have an x-axis as latitude and a y-axis containing the elevation attributed by the geolocated photons.Figure 7c shows the result of the DBSCAN algorithm, which successfully distinguished the photons returned from land, water surface, water column, and seafloor.By default, photons that have returned from the seafloor are apparent and need applying refraction correction to retrieve their actual depths 41 ; this is because there will be a change in the speed of light that occurs at the air-water interface due to the fact that the refractive index of air and seawater is different.In our research, we have implemented the refraction correction to the seafloor returned photons using the method suggested by Dandabathula et al. 29 , Parrish et al. 41 , and Guo et al.48 ; the formula used to perform the refraction correction for the photons that returned from the water column and the seafloor is shown in Eq. ( 1).
where Depth C is the refraction-corrected depth of the seafloor, Depth apparent is the apparent/default depth (without refraction correction) obtained by computing the difference between the elevation of the water surface (moving average of the surface waves) and the seafloor, and n 1 ≈ 1.00029 and n 2 ≈ 1.34116 are the refractive index of air and water bodies.The refracted corrected photons returned from the seafloor are shown in Fig. 7c.Photons returned from the water surface and water column were discarded during the computation of bathymetry as they do not have any role in estimating the seafloor depth.
The default vertical datum of ICESat-2 ATL03, i.e., heights above the WGS84 ellipsoid, was converted to orthometric heights by preferring the EGM2008 geoid model in the geoid height calculator available at https:// www.unavco.org.By extending this method to all the 396 beams over the study area, approximately 0.2 million ICESat-2 photons representing the seafloor's depth values and terrain elevations were collected as a part of data collection from the ICESat-2 ATL03 data products.Figure 8 shows the schematic representation of the methodology implemented in this research towards generating a DBEM for Adam's Bridge, which is partly modified from earlier researchers' works 29,35,41,[47][48][49]52,[86][87][88] .
Interpolation is essential in generating bathymetric surfaces from known depth points; the process enables estimating the depths in areas lacking direct measurements.The point database containing the accumulated depths and elevation values from the ICESat-2 photons must undergo an interpolation procedure to generate a continuous surface.Due to the along-track pattern of data collection by ICESat-2, interpolation techniques may not yield a continuous surface; for this, additional points representing seafloor depths are needed in the regions with gap areas.The additional depths were considered from available ENCs and charts (mentioned in Table 1).Similarly, to increase the point density over the extent of land (Rameshwaram and Talaimannar islands), additional elevation values (converted to orthoheights) were considered from FABDEM 58 .Figure 9 shows the distribution of seafloor depths accumulated from the processed ICESat-2 photons and ENCs/charts, along with the elevation values for Islands from FABDEM.A few points from the ENCs/charts were reserved for technical validation (considered as checkpoints) of the output DBEM (discussed in the subsequent section) and not included as a part of the above-said point database.

Validation of depths derived from ICESat-2 water penetrated photons
Towards assessing the accuracy of the depths obtained from the ICESat-2 water-penetrated photons, technical validation has been performed at two test sites that are falling in the study area of this research (shown in Fig. 10).High accurate reference data is needed to compare the depths obtained from ICESat-2 photons.However, due to the shallow nature of the water around the Adam's Bridge, no research vessel could attempt to negotiate the dangers associated with the area.To date, geological studies around the Adam's Bridge are meager except for a limited geomorphological observation made by the Marine Wing of the Geological Survey of India (GSI) at the southern part of Dhanushkodi Island 85 and a seismic survey by the National Institute of Oceanography, Council of Scientific and Industrial Research at the northern extent of Dhanuskhodi Island 68 .In particular, GSI, as a part of the geomorphological observation, conducted a detailed bathymetric survey for an area of about 12 × 2 km by engaging a mechanized wooden boat fitted with a portable echo sounder 85 .Bathymetric contours resulting from this survey were used as reference data for test site 1 (refer to Fig. 10a).Similarly, i-Boating: Marine Navigation offers marine charts through its application that offers offline nautical charts, inland river navigation charts, and lake contour maps for fishing, kayaking, yachting, and sailing (available as a subscription service at https:// fishi ng-app.gpsna utica lchar ts.com/).The navigation chart offers a warning contour at 4 m depth in shallow water near the extent of Adam's Bridge and Talaimannar (refer to Fig. 10b); this warning contour at 4 m depth prompts large draught boats not to navigate beyond the safety contour.This warning contour at a depth of 4 m is used   2) at these 25 locations using Root Mean Square Error (RMSE), a statistical quantifier of error.Table 2 shows the depth values derived from the ICESat-2 water-penetrated photons and their differences with respect to the reference data.Root Mean Square Error (RMSE) shows how far predictions fall from true values using Euclidean distance.The formula used to compute RMSE is shown in Eq. (3).
where (Depth) ICESat-2 is the set of depth values obtained from the ICESat-2 water-penetrated photons and (Depth) reference is the set of depths considered from the reference datasets.n (= 25) is the number of observations.RMSE resulting from these 25 validation points is 0.2 m.Earlier, researchers performed the accuracy assessment of the seafloor depths derived from the ICESat-2 photons for shallow water, and it is proven that the performance is as per the sensor specifications 52 , and generally, during clear water conditions, the accuracy in terms of RMSE is in the range of 0.2-0.8m 41,45,49 .These accuracies can be termed significantly high accurate   measurements from the currently available operational space-borne active sensor systems.Moreover, the accuracy assessment from our research, which yielded an RMSE of 0.2 m, is also in the similar lines of assessment done by earlier researchers, indicating significant potential for bathymetric mapping with ICESat-2 photons, especially for the shallow waters.

Generation of DBEM and its accuracy assessment
The point database containing depths and the elevations of the study area were interpolated using the Inverse Distance Weighted (IDW) method; this method is considered a highly adaptable estimation method as it is best to reconstruct natural surface recourses given dense and well-distributed points in the study area 93,94  grids represent higher resolution and more detail.Determination of the cell size of a DEM is a fundamental component during a surface generation process as it aims to produce the best representation of a terrain surface without introducing artifacts from interpolation.A standard method for determining the cell size of a rasterbased elevation model has been defined by Hu 95 and endorsed by Langridge et al. 96 is shown in Eq. ( 4) where s is the estimated raster cell size, and n is the number of points in the minimum area of density (A) within the extent of the point distribution.From the point distribution map obtained in our experiment, it is observed that for every 200 sq.m, at least 2 points representing depth/elevation exist.Thus, the output pixel's cell size was set to 10 m during the interpolation stage based on Eq. ( 4).By applying the processing scheme as mentioned above, we generated the high-resolution DBEM for the Adam's Bridge, in which ~ 0.2 million points are from the ICESat-2 photons have played a significant role in the surface generation along with the additional points from ENCs/charts and FABDEM (elevation points for Islands).During this process, a few values of sounding depths (considered from INHO chart number 3040) were reserved as checkpoints towards quality checking of the output DBEM, which were not used in the interpolation.

Figure 1 .
Figure 1.Location map of Adam's Bridge.(a) The map shows India and Sri Lanka by highlighting the extent of Adam's Bridge with a yellow box.(b) The map shows Rameshwaram Island (Indian side) and Mannar Island (Sri Lankan side) and highlights Adam's Bridge with a white box that constitutes the study area.(c) Satellite Image showing Adam's Bridge, which stretches from Dhanushkodi to Talaimannar.Source Map prepared using the satellite imagery through web mapping services (https:// tiles.maps.eox.at/ wms)* of the Sentinel-2 cloudless layer for 2021 by EOX ( and ).Maps were composed using QGIS Ver.3.30.0-'s-Hertogenbosch().*This work is licensed under a Creative Commons by Attribution (CC BY-NC-SA 4.0) license.

Figure 2 .
Figure 2. Three-dimensional perspective views of digital bathymetric elevation model generated using ICESat-2 water-penetrated photons for Adam's Bridge.(a) Perspective view of Adam's Bridge from the Gulf of Mannar as an observer position.(b) Perspective view of Adam's Bridge from Palk Strait as an observer position.The primary inference from both these perspectives is that in its entire form, Adam's Bridge is a submerged ridge with a submarine continuation of Dhanushkodi and Talaimannar Island.The exposed sand banks are only 0.02 percent of Adam's Bridge's total volume when the base is considered at 8 m depth.These perspectives were generated using ESRI's ArcScene Ver.10.8.1 software (https:// www.esri.com/).

Figure 3 .
Figure 3. Contour map of 2 m interval for the extent of Adam's Bridge derived from the digital bathymetric elevation model generated using ICESat-2 water-penetrated photons.The topographic surface on either side for 1.5 km of the Adam's Bridge's crest line is highly undulating with sudden depths.Towards the Gulf of Mannar side, from Dhanushkodi's tip, the distance from the crest of Adam's Bridge to the − 8 m contour line is 6.5 km for its initial 17 km.Towards reaching the Talaimannar side, the distance to the base fluctuates between 2.5 and 3.5 km.Towards the Palk Strait side, the distance from the crest to the base is always greater than 7 km.

Figure 4 .
Figure 4. Slope map for the extent of Adam's Bridge derived from the digital bathymetric elevation model generated using ICESat-2 water-penetrated photons.The transverse slope towards the Gulf of Mannar mainly varies between 0.05° and 0.35°.In contrast, towards the side of Palk Strait, the slope, on average, is never greater than ~ 0.2°.

Figure 5 .
Figure 5. Depiction of the narrow channels along Adam's Bridge's crest line.(a) Satellite image showing the exposed sand banks along the Adam's Bridge crest line.(b) The Elevation profile along Adam's Bridge's crest line shows narrow channels up to a depth of 3 m.(c) The three-dimensional perspective of Adam's Bridge shows the narrow channels enabling water exchange between the Gulf of Mannar and the Palk Strait.Source Map in (a) is composed in QGIS Ver.3.30.0-'s-Hertogenbosch(https:// qgis.org) using the satellite imagery through web mapping services ()* of the Sentinel-2 cloudless layer for 2021 by EOX ( and ).The perspective view in (c) was generated using ESRI's ArcScene Ver.10.8.1 software ().*This work is licensed under a Creative Commons by Attribution (CC BY-NC-SA 4.0) license.

Figure 6 .
Figure 6.ICESat-2's reference ground tracks and available beams over the extent of Adam's Bridge.(a) Seven reference ground tracks of ICESat-2 are available over the extent of Adam's Bridge.(b) Available beams (both strong and weak) over the extent of Adam's Bridge.Source Maps are composed in QGIS Ver.3.30.0-'s-Hertogenbosch(https:// qgis.org) using the satellite imagery through web mapping services ()* of the Sentinel-2 cloudless layer for 2021 by EOX ( and ).*This work is licensed under a Creative Commons by Attribution (CC BY-NC-SA 4.0) license. https://doi.org/10.1038/s41598-024-65908-2

Figure 7 .
Figure 7. Detection of seafloor from ICESat-2 water-penetrated photons.(a) A subset of the ICESat-2 beam acquired over the extent of Adam's Bridge.(b) The plot shows the geolocated photons, including noise, that are recorded in the Level-2A ATL03 data product.(c) Classification of ICESat-2 photons based on the DBSCAN algorithm; the result consists of photons returned from land, water surface, water column, and seaflooor.The ray tracing mechanism of ICESat-2 photons in the water column will be impacted due to the difference in air and water refraction index as it causes apparent depths; the refraction correction procedure significantly improves the depth estimation accuracy.Source Map in (a) is composed in QGIS Ver.3.30.0-'s-Hertogenbosch(https:// qgis.org) using the satellite imagery through web mapping services ()* of the Sentinel-2 cloudless layer for 2021 by EOX ( and ).*This work is licensed under a Creative Commons by Attribution (CC BY-NC-SA 4.0) license.

Figure 9 .
Figure 9. Distribution of points representing the depths accrued from seafloor returned, and refraction corrected ICESat-2 photons, sounding depths from ENCs/charts, and land elevation values from FABDEM.

Figure 10 .
Figure 10.Test sites to evaluate the bathymetric performance of ICESat-2 geolocated photons.(a) Test site 1 shows the extent of southern Dhanushkodi with GSI-surveyed bathymetric contours overlaid with the ICESat-2 beams.(b) Test site 2 shows the extent of Adam's Bridge and Talaimannar with i-Boating: Marine Navigationbased 4 m shallow water warning line overlaid with the ICESat-2 beams.Source Maps are composed in QGIS Ver.3.30.0-'s-Hertogenbosch(https:// qgis.org) using the satellite imagery through web mapping services ()* of the Sentinel-2 cloudless layer for 2021 by EOX ( and ).*This work is licensed under a Creative Commons by Attribution (CC BY-NC-SA 4.0) license.

Figure 11 .
Figure 11.Accuracy assessment of digital bathymetric elevation model (DBEM) generated with collection of depths accrued from ICESat-2 geolocated photons and ENCs.(a) DBEM showing the test site (marked in magenta) containing locations of validation points (b) Chart showing sounding depths and locations of validation points.

Table 1 .
Details of the data sources used to generate a digital bathymetric elevation model for the Adam's Bridge.Elevation values from FABDEM were used to increase the density of points in the study area for the extent of the land part (Rameshwaram and Talaimannar Islands) Preferred beams: Higher preference for the strong beams but cannot omit the photons in the weak beams when there is a need to accumulate denser points to generate a bathymetric surface Post-processing methods: Refraction correction applied to the photons returning from the seafloor Total number of accrued depths from seafloor returned photons: 0.2 million K d (490) layer from Sentinel-3 A/B mission Source: https:// senti nels.coper nicus.eu Usage: For the dates acquired by the ICESat-2, within + / − 24 h, corresponding Level-2 Ocean and Land Colour Instrument (OLCI) data products from the Sentinel-3 A/B mission were used to retrieve K d (490).K d (490) data was used to assess the turbid load in the study area.Only those acquisitions of ICESat-2 were considered while K d (490) < 0.12 m −1 , i.e., clear water conditions Sounding depths from Electronic Navigational Charts (ENCs) and charts Source information of ENCs: https:// hydro bharat.gov.in, http:// www.nara.ac.lk, and https:// iho.int Usage: The sounding depth values were used to increase the density of depth values in the study area Procedure adopted: Sounding depths were digitized from ENCs and chart datasets issued by hydrographic offices Forest And Buildings removed Copernicus DEM (FABDEM) Source: https:// data.bris.ac.uk Usage: Sounding depths from chart number 3040.(Surveyed by Indian Naval Hydrographic Office) Available at an extent containing the junction of Dhanushkodi Island and Adam's Bridge Source: https:// hydro bharat.gov.in/ produ cts/ online-catal ogue/ Usage: Used as reference data at 10 locations to assess the accuracy of DBEM generated in this research Vol.:(0123456789) Scientific Reports | (2024) 14:14896 | https://doi.org/10.1038/s41598-024-65908-2

Table 2 .
Values of depths from water penetrated ICESat-2 photons and reference datasets.

Table 3 .
Comparison of depth values retrieved from Indian National Hydrographic Office chart and ICESat-2 photons based digital bathymetric elevation model.